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We investigate the Ruderman-Kittel-Kasuya-Yosida oscillations of the itinerant carrier spin den- 
sity in a system where those oscillations appear only due to a finite distribution of a localized spin. 
The system represents a half-infinite one-dimensional quantum wire with a magnetic impurity lo- 
cated at its edge. In contrast to the conventional model of a point-like exchange interaction the 
itinerant carrier spin density oscillations in this system exist. We analytically demonstrate that 
when the radius of the exchange interaction is less than the wave length of the itinerant carriers 
living on the Fermi surface, the long range behavior of the oscillations is identical to the one taking 
place in the zero radius limit of the same exchange interaction but for an infinite one-dimensional 
quantum wire where, in comparison with the original half-infinite system, the mass of the itinerant 
carriers is strongly modified by the exchange interaction radius. On the basis of our analysis we 
make a suggestion on directionality of surface Ruderman-Kittel-Kasuya-Yosida interaction shown 
in recent experiments: we believe that in general the anisotropy of the Ruderman-Kittel-Kasuya- 
Yosida interaction could result not only from the anisotropy of the Fermi surface of itinerant carriers 
but also from the anisotropy of the spin carrying atomic orbitals of magnetic impurities. 

PACS numbers: 75.20.Hr, 75.75.-c, 71.10.-w 



I. INTRODUCTION 

Systems with interacting localized magnetic moments 
play an important role in modern spintronicsA. Among 
those systems the ones where the interaction between lo- 
calized spins is due to the itinerant carriers are very at- 
tractive for semiconductor spintronics because the mag- 
netic properties of the systems can be manipulated by 
the itinerant carrier density. The latter in semiconduc- 
tors can be easily varied using external gate voltage. 

This, so-called itinerant magnetism, is provided by 
the exchange interaction introduced by Ruderman and 
Kittel^ in the context of the hyperfine interaction be- 
tween a localized nuclear spin and conduction electrons 
and later by Kasuya and Yosida&A in the context of the 
s — d interaction between conduction electrons and mag- 
netic ions of transition metals. 

The essence of the Ruderman-Kittel-Kasuya-Yosida 
(RKKY) itinerant magnetism is that the exchange in- 
teraction between localized spins and paramagnetic itin- 
erant carrier system induces oscillations of the itinerant 
carrier spin density. These oscillations propagating in 
an originally non-magnetic system couple localized spins 
and create a magnetic structure. This type of interac- 
tion becomes relevant when the direct spin-spin coupling 
is weak, i.e., in diluted systems. For example ferromag- 
netism in bulk diluted magnetic semiconductors^ is of the 
RKKY type. It was investigated in Gai-^Mn^As using 
high-resolution scanning tunneling microscope^. 

The RKKY effect has also been widely studied in meso- 
scopic systems both without spin-orbit interactions^— 
and with spin-orbit interactions^—. In the latter case 
spin-orbit interactions make the RKKY spin density os- 
cillations anisotropic. Impact of electron-electron inter- 
actions on the RKKY range function in one-dimensional 
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FIG. 1: (Color online) A half-infinite ID QW with a magnetic 
ion placed at its edge taken as a reference point. The ionic 
spin is S and ro characterizes its spatial distribution. The 
electronic density is zero at the edge, x = 0. Thus for ro — 
one expects no RKKY oscillations in the system. However, for 
x > the electronic density exhibits the Friedel oscillations 
converging to the electronic density of the corresponding in- 
finite ID QW, no. Therefore for ro 7^ the exchange interac- 
tion takes place and is responsible for the RKKY oscillations 
of the electron spin density. 



quantum wires (ID QW) was considered in Refs. 
It was shown that the range function is different from the 
one in non-interacting ID QWs. It turns out that due to 
electron-electron interactions the range function decays 
slower than in the non-interacting case. 

However, in spite of considerable amount of scientific 
research on the RKKY interaction a little attention has 
been paid to the fact that the exchange interaction re- 
sponsible for the RKKY interaction is actually non-local 
in space but has a finite range ro given by the spatial 
spin distribution around the localization center. In the- 
oretical models of the RKKY interaction it is almost al- 
ways assumed that ro = and the exchange interaction 
has a point-like character which is modeled by the Dirac 
delta-function. Of course, in many cases it is a very good 
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approximation when the wave function of a magnetic im- 
purity is well localized which means that r is less than 
all other relevant lengths present in the system. How- 
ever, in semiconductors there is one important length, 
namely the wave length of the itinerant carriers living 
on the Fermi surface, If, which can be varied in a wide 
range. This can be done, e.g., by changing the carrier 
density with the help of a gate voltage. Therefore, one 
can reach a situation where If < rg. 

What is more interesting is that even in the case when 
If > ?'o h can happen that in a given system the RKKY 
spin density oscillations appear only due to a finite value 
of r and cannot appear if r = 0. Such a situation 
takes place, for example, when magnetic impurities are 
deposited on a surface of a paramagnetic bulk sample. 
The wave function of the itinerant carriers in the sam- 
ple is equal to zero on its surface and therefore for a 
point-like model of the exchange interaction, ro — 0, the 
RKKY spin density oscillations do not appear and can- 
not couple the spins of the magnetic impurities on the 
surface. Therefore, no magnetic ordering will take place. 
However, in reality ro =/= and the exchange interaction 
will definitely produce the spin density oscillations with 
the subsequent magnetic structure of the impurity spins 
on the surface of the sample. Thus, in this system a finite 
value of ro creates itinerant surface magnetism. 

Some efforts to consider a finite range of the exchange 
interactions were made in the research devoted to ferro- 
magnetism in Heusler alloys^—. An essential drawback 
of the models of the exchange interactions in those at- 
tempts was that in the limit ro — > the models did not 
reproduce point-like exchange interactions. However, it 
is desirable to reproduce the Dirac delta-function in the 
limit ro — > because it is often used in many theoretical 
models and therefore one would like to make a compari- 
son with traditional theories in the limit ro — > 0. 

A finite range exchange interaction which in the limit 
ro —> gives the traditional Dirac delta-function model 
has been used in Ref. |23to calculate the RKKY exchange 
integrals and in Ref. [23| to study the RKKY spin density 
oscillations in a mesoscopic ring. In particular, in Ref. [22j 
it has been demonstrated that a finite value of ro removes 
an unphysical divergency of the RKKY integral at zero 
distance while in Ref. l23| it has been shown that any 
finite value of ro always becomes important if the number 
of the electrons in the ring is large enough. However, 
the research in Ref. [23| represented mainly a numerical 
experiment and, as a result, did not clarify the physical 
role played by the localization radius of the impurity spin. 
In addition, the RKKY spin density oscillations in that 
study were also present in the limiting case ro — > and 
therefore a finite value of ro did not play any key role in 
the formation of the RKKY spin density oscillations. 

In the present work we use the model of Ref. [H for a 
finite range exchange interaction in a system where the 
RKKY spin density oscillations appear only if ro 7^ 0. 
The system represents a half-infinite ID QW with a mag- 
netic impurity on its edge. Since the wave function of the 



itinerant carriers (electrons to be definite) is zero on the 
edge of the wire, the RKKY spin density oscillations ex- 
ist only if ro is finite. To understand the physical role 
played by the impurity spin localization radius in the 
RKKY spin density oscillations we solve the problem an- 
alytically in the case when lp > tq and calculate the 
coordinate dependence of the RKKY oscillations in this 
system. We demonstrate that at large distances this de- 
pendence is identical to the one of the RKKY oscillations 
which one would obtain in the limiting case r = of our 
model but for an infinite ID QW with electrons whose 
mass is strongly modified by the localization radius of the 
impurity spin. 

The paper is organized as follows. In Section [H] we 
mathematically formulate the problem and solve it in 
Section [Til] using the Feynman diagram approach. In 
connection with recent experiments we make a sugges- 
tion in Section IIVI about directionality of surface RKKY 
interaction. Conclusions are given in Section [Vj 



II. FORMULATION OF THE PROBLEM 

As mentioned in the introduction the system represents 
(see Fig. [T|) a half- infinite ID QW with, e.g., a magnetic 
ion on its edge. It could be also another object, natural or 
artificial (e.g., a quantum dot with odd electron number) 
with non- vanishing total spin S. The electrons in a half- 
infinite ID QW have eigenenergies 
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where m is the effective mass of the conduction electrons 
and q x is not the electron momentum because the trans- 
lational invariance is broken. The quantum number q x 
takes all real values except zero, q x ^ 0. The eigenstates 
\q x cr), where a is the spin quantum number, in the coor- 
dinate representation are 



[xa \q x a) = 



&<y'<?-jl sin(ig a; x) x ^ 
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where L is the size of the system. To consider a half- 
infinite system we take the limit L — > 00 in the sub- 
sequent calculations. Since the momentum operator is 
p x = —ihd/dx, one can see that (xa'\q x a) is not the 
eigenfunction of p x and thus, as mentioned above, the 
quantum number q x is not the electron momentum. 

The electronic density vanishes on the edge of the wire. 
Far from the edge it has the value which it would have 
in an infinite ID QW with the Fermi momentum qp (in 
the infinite case it would be real momentum), 
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In the region close to the edge the electronic density 
shows oscillations known as the Friedel oscillations. 
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The half-infinite quantum wire can be modeled by the 
following potential: 



u(x) 
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(3) 



The second quantized Hamiltonian of the half-infinite 
ID QW without the magnetic impurity thus reads in the 
coordinate basis as 
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where ^J(x), ipv{ x ) are the electronic creation and anni- 
hilation field operators of the half-infinite ID QW. 

When the magnetic impurity is deposited on the edge 
of the half-infinite ID QW, the electronic spin density 
interacts through the exchange interaction with the im- 
purity spin. The corresponding second quantized Hamil- 
tonian is 



H„=Y,S i (tT\a i \o') J dxJ(x)^l(x)^(x), 



(5) 



where a is the vector of the Pauli matrices and for J(x) 
we use the model of Ref. n29, 
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The total Hamiltonian describing the RKKY electronic 
spin density oscillations is given by the sum 



(7) 



To get the RKKY oscillations and understand the 
physical role of the finite impurity spin distribution ro 
on those oscillations one may consider H ex as a pertur- 
bation. In the next section we will implement this pertur- 
bation theory on the language of the Feynman diagram 
expansion. 



III. DIAGRAMMATIC SOLUTION 

To obtain the electronic spin density we first have 
to find the imaginary-time (or Matsubara) one-particle 
Green's function, 

g^(xr\x'T') = (TMx,r)i>t t (x',r')}, (8) 

where the angular brackets stand for the thermal average 
and the electronic field operators are in the imaginary- 
time Heisenberg representation, 
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(9) 
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FIG. 2: (Color online) The first order Feynman diagram con- 
tributing to the electronic spin density oscillations. Unlike 
the traditional way to draw the propagators by straight lines 
the incoming and outgoing propagators are drawn using wavy 
lines to stress that they take into account the Friedel oscil- 
lations of the electronic density which are present in a half- 
infinite ID QW without the magnetic impurity. 



where (i is the chemical potential. 

The exact electronic spin density is obtained from the 
imaginary-time Green's function as 



ai xact (x) = -Ti[a i g(xr\xr + 0)}, 



(10) 



where the trace is taken in the spin space. 

The RKKY oscillations of the electronic spin density 
are contained in the contribution to Q aa '(xT [x't 1 ) of the 
first order in the perturbation Hamiltonian H ex . The 
corresponding Feynman diagram is shown in Fig. [5] The 
wavy lines in this diagram are the electronic propaga- 
tors g ( °}, (xt\x't') of the half-infinite ID QW without the 
magnetic impurity. They include the Friedel oscillations 
of the electronic density and can be found using a trick 
suggested in Ref. IU The trick is based on the image 
charge technique- -. The result is 

g { 3(xr\x'r') = g ( £{xT\x'T') - gW(xr\ - x't') (11) 

for x and x' 0, 

g^,(xT\x'r') = (12) 

for x < or x' < 0. In Eq. (Till ^°j,(arr|arV) is the 
imaginary-time Green's function of the corresponding in- 
finite ID QW without the magnetic impurity. 

Using the rules- - for the analytic reading of Feyn- 
man diagrams we obtain from the diagram shown in 
Fig. [5]the first order contribution to g aa '(xT\x r T'), i.e., 

g^}i{xT\x' t') and calculate the electronic spin density 



a l (x) = -Tr[a l g {1 \xT\xT + 0)]. 



(13) 



It is easily verified that the result is given by the following 
expression: 

. , . 16m J S 1 

a 2 

x J dq x J d^sinQ^a^sinQ^xJ x (14) 
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where x ^ 0, n„ x is the Fermion occupation number, 
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with /? = 1/ksT being the inverse temperature and the 
matrix Q 9x9 ' is 



e 4" 



(16) 



Therefore the dependence of the electronic spin density 
on the impurity spin distribution comes through the ma- 
trix Q fe9i . 

One can see from Eq. (TTB1) that when tq — 0, the 
matrix Qq x q' vanishes and no RKKY spin density oscil- 
lations appear in the system. This is in agreement with 
our earlier qualitative discussions. 

Below we only consider the most interesting case of 
zero temperature T — 0. Using Eq. (|16p we obtain from 
Eq. (|14[) in the zero temperature limit the following ex- 
pression for the electronic spin density: 

u , 8mJS i 

I dq x I dq x sm(^^q x \sin(^^q x jx ^ 
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^ffe-C) 2 
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ql - & 



where q x and q' x are dimensionless integration variables. 
The integral 
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where lp — fr/qF, can be calculated analytically in the 
case ro/lp <C 1. To do this we perform the expansion 



^(^fe^) 2 .,^^) 2 ^) 2 ,,^ (19) 

Let us denote through J(x) the integral I(x) approxi- 
mated using this expansion, 



J(x) = x 
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The analytical expression for the integral J(x) is 
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FIG. 3: (Color online) Comparison of the oscillation integrals 
I(x) and J(x) for to/If = 0.4. The maximal deviation of the 
approximate integral J(x) from the exact integral I(x) takes 
place at short distances and is around 14%. 
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FIG. 4: (Color online) Comparison of the oscillation integrals 
I(x) and J(x) for to/If = 5. The approximate integral J(x) 
greatly overestimates the exact integral I(x). Moreover, the 
short range behavior of J(x) is even incorrect qualitatively. 
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To demonstrate that the integral J(x) is really a good 
approximation to the exact integral I{x) we calculate the 
integral I(x) numerically and show two situations in Figs. 
[3] and |4] for the two cases ro/lp < 1 and ro/lp > 1, 
respectively. Even when r^/lp = 0.4, i.e., is not too much 
smaller than 1, the integral J(x) gives a good precision 
and for smaller ro/lp the difference between J(x) and 
I(x) decreases rapidly. On the other side, as one can 
see from Fig. HI in the case r /lp = 5 the integral J(x) 
gives values much larger than the exact integral I{x) and 
at small distances it even predicts a wrong qualitative 
coordinate dependence. 

In the case when ro/lp < 1 we therefore obtain the 
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RKKY oscillations of the electronic spin density, 



a i (x)=JS l 



1 2m 
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(22) 



In diluted magnetic systems one is interested in the 
long range behavior of the RKKY oscillations. From Eq. 
(1221) it follows that the long range behavior of the RKKY 
oscillations is 



a l (x) = JS l 



' 2m cos(2x/l} 
^¥ 2x/l F 



(23) 



Comparison of the long range behavior of the RKKY 
spin density oscillations, Eq. (|23[) . with the one in a 
point-like Dirac's delta-function model of the exchange 
interaction in an infinite ID QW— allows us to make 
the following conclusion. The RKKY oscillations induced 
by finite range exchange interactions, ro ^ 0, in a half- 
infinite ID QW are identical to the ones which would take 
place in the limit ro = of our model but in an infinite ID 
QW with conduction electrons having a modified mass 



(24) 



This modification is quite strong because in semiconduc- 
tors the quantity r^/lp may be tuned in a wide range 
and it can be made much less than 1. 

Eq. (|M|) clarifies the physical role of the finite mag- 
netic impurity spin distribution in the formation of the 
RKKY spin density oscillations. One can still use a point- 
like Dirac's delta-function model of the exchange inter- 
action between the magnetic impurity spin and conduc- 
tion electrons. However, the electron mass m must be 
replaced by the modified mass m*. 



IV. A SUGGESTION ON DIRECTIONALITY 
OF SURFACE RKKY INTERACTION 

In recent experiments^ surface RKKY interaction was 
studied using cobalt adatoms on platinum (111). A 
strong directional dependence of the RKKY spin density 
oscillations was observed. The anisotropic behavior was 
attributed only to the anisotropy of the Fermi surface of 
platinum. 

In connection with our analysis we would like to note 
that since the wave function of the itinerant carriers in 
platinum is equal to zero on its surface, the RKKY spin 
density oscillations in the experiments of Ref. [28| are 
definitely only due to a finite spin distribution of cobalt 
adatoms. As we have shown above for a simple one- 
dimensional model, the RKKY spin density oscillations 
existing only due to a finite value of ro are created by 
electrons whose mass is strongly modified by the radius of 
the exchange interaction rp. Therefore, in addition to the 



Fermi surface anisotropy the anisotropic behavior of the 
RKKY spin density oscillations could also be produced 
by the fact that the itinerant carrier mass for different 
directions suffers different modification. 

It is not difficult to generalize our model, Eq. (|6|), 
to take into account the anisotropy of the spin carrying 
atomic orbitals of a magnetic impurity. The simplest 
possibility for the case of magnetic adatoms on a two- 
dimensional surface is 
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(25) 



where r$ ^ Tq. Other more adequate models taking into 
account the real angular distribution of the magnetic im- 
purity spin are also possible. 

Since in general the symmetries of the Fermi surface of 
itinerant carries and the spin carrying atomic orbitals of 
a magnetic impurity are different, it is experimentally 
challenging to distinguish these two different kinds of 
symmetries in the RKKY spin density oscillations. The 
experimental observation of traces of the symmetry of 
the spin carrying atomic orbitals of a magnetic impurity 
could prove that a finite spin distribution of the impu- 
rity spin plays a significant role in the formation of the 
RKKY interaction. 



V. SUMMARY 

We have studied the Ruderman-Kittel-Kasuya-Yosida 
(RKKY) spin density oscillations in a system where these 
oscillations appear only due to a finite spin distribution of 
a magnetic impurity. It has been shown that the physical 
role played by the finite spin distribution is to modify the 
mass of the itinerant carriers in the system. 

Making use of this result we have made a suggestion on 
the anisotropic behavior of surface RKKY interaction ob- 
served in recent experiments which were explained only 
from the point of view of the anisotropy of the Fermi sur- 
face of itinerant carriers. We assume that the anisotropy 
of the spin carrying atomic orbitals of a magnetic im- 
purity could result in an anisotropic itinerant carrier 
mass modification and thus it could also contribute to 
an anisotropic behavior of surface RKKY interaction but 
with the impurity atomic orbital symmetry which is in 
general different from the symmetry of the Fermi surface 
of the itinerant carriers. 
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